Quantifying non-Markovianity of continuous variable Gaussian dynamical maps 
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We introduce a non-Markovianity measure for continuous variable open quantum systems based 
on the idea put forward in H.-P. Breuer et al. Phys. Rev. Lett. 103, 210401 (2009), i.e., by quanti- 
fying tiie flow of information from the environment back to the open system. Instead of the trace 
distance we use here the fidelity to assess distinguishability of quantum states. We employ our 
measure to evaluate non-Markovianity of two paradigmatic Gaussian channels: the purely damping 
channel and the quantum Brownian motion channel with Ohmic environment. We consider different 
classes of Gaussian states and look for pairs of states maximizing the backflow of information. For 
coherent states we find simple analytical solutions, whereas for squeezed states we provide both 
exact numerical and approximate analytical solutions in the weak coupling limit. 
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I. INTRODUCTION 

Physical systems are never perfectly isolated from their 
environment. Especially in quantum mechanics, when 
they are exploited to perform quantum computation or 
communication protocols, the interaction of the systems 
of interest with the environment should be considered in 
the derivation of the dynamical equations. Effects of this 
interaction, e.g. decoherence and disentanglement, can 
indeed endanger the accomplishment of any task based 
on quantum features like coherence or entanglement. 

In order to take into account the presence of the envi- 
ronment and its influence on the dynamics of the system, 
within the theory of open quantum systems P, I3 a va- 
riety of techniques have been developed to describe the 
evolution of the system of interest, e.g. by quantum mas- 
ter equations. The functional form of any master equa- 
tion depends both on the system and environment, and 
on the specific features and strength of the interaction. In 
the literature the dynamics of open quantum systems are 
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often described using master equations in the so-called 
Lindblad form Q . Profitable applications of this class of 
dynamical equations are present in many fields of physics, 
and systems whose dynamics is described by equations in 
the Lindblad form are generally called Markovian. As it 
will be explained in more detail in Sec. II, the dynami- 
cal maps associated with Lindblad master equations are 
divisible, implying that during the evolution any pair of 
initially different states becomes less and less distinguish- 
able. This phenomenon is interpreted as an irreversible 
loss of information which flows from the system to the 
environment, and it is considered to be the key feature 
of Markovianity [3] • 

In practice, however, Lindblad master equations are 
derived under a series of approximations. The exact dy- 
namics of any physical system is generally described by 
other classes of master equations. Recently, a lot of ef- 
fort has been devoted to provide a formal definition of 
non-Markovianity in open quantum systems, e.g. to cap- 
ture physical features such as the re-coherence due to 
reservoir memory effects [1, 0. These efforts d, 0, II] 
also lead to computable measures for the degree of non- 
Markovianity. In this paper, we focus on the definition 
given in Ref. Q where non-Markovianity is defined in 



2 



terms of the information flow between the open system 
and its environment. 

Besides its own importance from a purely theoretical 
point of view, the concept of non-Markovianity and its 
quantification may also find practical applications. One 
may ask indeed, whether non-Markovianity can be con- 
sidered as a resource to improve quantum technologies. 
More specifically, assuming that the density of modes of 
the reservoir may be engineered in a controlled way to 
induce non-Markovian behavior, can this be used to im- 
prove existing quantum protocols? The first affirmative 
answers come from quantum metrology and quantum key 
distributon. In Ref. Q the authors investigate the prob- 
lem of parameter estimation when the quantum chan- 
nel is non-Markovian according to the definition given 
in Ref. Q. They find that, for some non-Markovian 
reservoirs, the estimation can be improved compared to 
the Markovian case. The other example is reported in 
Ref. ^lOi] where it has been proven that quantum key dis- 
tribution protocols in non-Markovian channels provide 
alternative ways of protecting the communication which 
cannot be implemented in usual Markovian channels. 

Even if the definition introduced in Ref. [4| is indepen- 
dent of the nature of the physical system, it has been 
applied so far to the discrete variable case only. In this 
paper we extend the analysis to continuous variable (CV) 
systems ll| which, in quantum information and commu- 
nication, represent a valid, and sometimes better, alter- 
native to discrete variable systems. Our aim is to in- 
troduce and study a computable measure for the degree 
of non-Markovianity in continuous variable systems fo- 
cussing on some relevant examples of Gaussian preserving 
maps. Moreover we also consider the possibility of eval- 
uating the map only for subsets of Gaussian states (e.g. 
coherent states and squeezed states) with the main intent 
to provide a characterization of the map for protocols 
relying only on those specific classes of states. This ap- 
proach paves the way to a definition of non-Markovianity 
as a resource in quantum information theory. 

The paper is organized as follows: in Sec. II we re- 
view the non-Markovianity measure we use and extend 
the definition to the realm of continuous variable Gaus- 
sian states. In Sec. Ill we focus on a phenomenological 
master equation describing a damping channel and eval- 
uate its non-Markovianity, whereas in Sec. IV we address 
the same issue for quantum Brownian motion in the weak 
coupling limit. Finally in Sec. V we discuss the results 
and close the paper with some concluding remarks. 



II. QUANTIFYING NON-MARKOVIANITY IN 
CONTINUOUS VARIABLE SYSTEMS 



The measure for the degree of non-Markovianity of a 
quantum process introduced in is based on the distin- 
guishability of two different initial quantum states pi and 
P2 under the action of the open system dynamical map 
<I>t associated to the process. The distinguishability is 



qualified and quantified in through the introduction 
of a proper distance measure between quantum states, 
the trace distance defined as I?(pi,p2) = Tr|pi — P2I/2. 
The trace distance satisfies a contractivity property un- 
der the action of any completely positive (CPT) map <I> 



(1) 



Given a dynamical map '^t.ta^ this is called divisible if the 
evolution up to a time t can be written as a completely 
positive evolution from the initial time to an inter- 
mediate time T, and another completely positive evolu- 
tion from the intermediate time to the final time t, i.e. 
*i't,to = ^t,T • ^T-jn, for any to < t < t. Lindblad dynami- 
cal semigroups describe divisible processes. It follows 
that under such dynamics the trace distance is always 
monotonic, i.e. 'D(^pi{t), P2{t)^> 'D(^pi{t), p2{t)^ for any 
pair pi(0), /52(0) of initial states and for any to < t < t. 
The monotonic decrease however holds also for classes 
of divisible maps which are not of a Lindblad form (ex- 
amples are given in [3]). On the other hand, the most 
dynamical evolutions violate both the divisibility condi- 
tion and the monotonicity property of the trace distance. 
When monotonicity is not satisfied it means that there 
are intervals of time for which the states become more 
distinguishable compared to previous instants. This fea- 
ture is interpreted as a fiow of information from the en- 
vironment back to the system, a striking property which 
characterizes a non-Markovian evolution. The measure 
of the degree of non-Markovianity is then defined as 



-V{piit),p2{t))dt, 



(2) 



where V indicates the time derivative and the maximiza- 
tion is taken over all the possible pairs of initial states. 

So far the quantity in ^ has been evaluated and an- 
alyzed in some details for discrete variables quantum 
maps, e.g. a one qubit channel j3, In this paper 

we extend it to continuous variable systems, focusing to 
single-mode systems. The extension involves two main 
issues requiring specific attention. The first comes from 
the fact that the Hilbert space for continuous variable 
systems is infinite dimensional, and therefore it is not 
possible to characterize all the states with a finite number 
of parameters as in the qubit case (e.g. Bloch sphere rep- 
resentation) . The second issue is related to the lack of an 
analytic expression for the trace distance or other, equiv- 
alent, distance measures for a generic CV state. On the 
other hand, these issues may be solved upon restricting 
the analysis to Gaussian states, and Gaussian preserving 
channels 15|. In fact, Gaussian states can be uniquely 
characterized by a finite number of parameters. More- 
over, since analytic expressions for the trace distance are 
lacking, alternative distinguishability signatures may be 
employed within the same spirit. One possible choice is 
to use the fidelity 



-^(Pi, P2) = TrA/ 



(3) 
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which is related to a proper distance measure, the Bures 

distance 'Djr{pi,p2) = \J'^~ '^\/ ^{pii P'^- Remarkably, 
the fidelity, and, thus, also the Bures distance, are mono- 
tonic under the action of any CPT map, making it a good 
candidate for taking the role of the trace distance in the 
definition of the non-Markovian measure. 

An alternative choice we could consider is the relative 
entropy S{pi\\p2) — Tr[pilog/9i — pi log/52] whose ex- 
pression for Gaussian states is known [ij]. Despite the 
fact that also the relative entropy possesses a contractiv- 
ity property under CPT maps, it is not a proper distance 
measure (e.g. it lacks symmetry property), and also it is 
not boimded. For these reasons we base our study on the 
fidelity. 

The mostgeneral single-mode Gaussian state can be 
written as [l5[ 

p'' ^ D{P)S{0^,^{N)S\OD\P) , 

where S{C} = exp[i(eat2 - S^*a^)] and D{P) = exp[^at - 
(3* a)] are the squeezing operator and the displace- 
ment operator, respectively, and vth{N) = [NY "/(I + 
NY ""''^ is a thermal equilibrium state with N average 
number of quanta, a being the annihilation operator. 
Upon introducing the vector operator = (i?i,i?2) = 
where q = {a + a'^)/^/2 and p = {a^ — a)/{i\/2) are 
the so-called quadrature operators, we can fully charac- 
terize p'~^ by means of the first moment vector 

= {R^) = \/2(3fie[/3],5m[/3]) , 

where {A) = Tr[Ap], and of the 2x2 covariance matrix 
(CM) cr, with elements 



[T]hk = ^{RhRk + RkRh) - {Rh){Rk) = 1, 2 . 

The expression of the fidelity for a generic pair of Gaus- 
sian states pf, P2 can be given in a closed analytical form 
[TgI and, as expected, depends only on the values of the 
vectors JC^, and covariance matrices <Tj of the states in- 
volved 



^/AT5-VS 



where 



A = 4Dct(cri + 0-2), 
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Det((Ti) 



Det(cr2) 



, (4) 

(5) 
(6) 



The fidelity in Eq. Q is a function of all the parameters 
that characterize the pair of Gaussian states: two com- 
plex displacement amplitudes /3i = |/3i|e*^', two complex 
squeezing amplitudes = r^e**^' and two real thermal 
parameters Ni. In order to simplify the notation used in 



the following we introduce here a set of collective argu- 
ments 



Pn = {Ni,N2}, 
Ps = {ri,r2,(l)i,(j)2} , 

Pc = {|/3l|, 1/321,^1,^2}: 



(7) 



and denote, e.g. by J-'({Pc, Pn}, i) the fidelity between 
two mixed coherent states at time t of the evolution. The 
full set of parameters is denoted by the symbol P. 

A non-Markovianity measure may be obtained by in- 
tegrating the time derivative of the fidelity T{P,t) over 
the intervals in which it decreases. If the class of initial 
states is characterized by the set of parameters P the 
measure may be written as 



A/p = max 
p 



dt 



T{P,t)dt 



(8) 



where T indicates the time derivative and the maximiza- 
tion is taken over the set of parameters P. 

The measure A/p in ([8]) is obtained by maximization 
over the class of Gaussian states, a procedure which 
makes it particularly suitable to asses Gaussian preserv- 
ing channels. On the other hand, when applied to a 
generic channel, it cannot be considered as a global prop- 
erty. This is not a crucial issue for practical applications 
for at least two reasons. On the one hand, this choice al- 
lows to actually calculate and compare the degree of non- 
Markovianity for continuous variable channels, a task 
that would not be feasible for non-Gaussian states. On 
the other hand, it should be noticed that while in prin- 
ciple there are techniques to prepare any kind of single 
qubit states, the same does not apply to continuous vari- 
able systems. As a matter of fact the class of Gaussian 
states plays a crucial role in quantum information pro- 
cessing, since they can be characterized theoretically in a 
convenient way, and they can also be generated and ma- 
nipulated experimentally in a variety of physical systems, 
ranging from light fields to atomic ensembles. Therefore, 
our aim in the rest of the paper will be that of character- 
izing the Gaussian degree of non-Markovianity for some 
relevant Gaussian preserving channels. 

Once we restrict the investigation to Gaussian states 
we still have to face the problem of the maximization 
procedure, which may be challenging from the numerical 
point of view, since the domains of some of the involved 
parameters are unbounded. One way to deal with this is- 
sue is to focus on subclasses of Gaussian states, e.g. pure 
coherent states or squeezed states, and therefore reduce 
the number of parameters involved. On the other hand, 
it is also possible to bound the domain of definition of 
the parameters, invoking the same line of reasoning used 
previously: experimental accessibility. For example, in 
practice it is not possible to obtain an arbitrary amount 
of squeezing 17[. The consequence is then a limitation 



of the domain of definition and therefore a faster con- 
vergence of numerical maximization algorithms. In the 
cases we examine in the next sections however we will see 
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that it is not always needed to bound the domain of defi- 
nition of the parameters. This is because the maximizing 
pair of states depends on the strength of the interaction, 
i.e., the couphng constant, and for weakly coupled sys- 
tems, experimentally accessible values for the squeezing 
and displacement may characterize the maximizing pair. 

In the following Sections we will assume that the 
maximum is achieved for pure states, i.e., we assume 
iVi = N2 — and perform the maximization over the 
other parameters. This assumption may be proved for 
the case of coherent thermal states in the weak coupling 
regime, whereas we conjecture its validity for the other 
classes of states. 

In the next two sections we introduce two different ex- 
amples of master equations used to describe the dynamics 
of continuous variable systems and we study their non- 
Mar kovian Gaussian properties. 



the derivative of the fidelity with respect to time. Be- 
cause the evolution of any Gaussian state depends on 
time through the function x{t) only, we have 



dJ- dF dx 2 dF 
dt dx dt dx 



(13) 



If 7(i) = 7o > and therefore XM{t) — 2a7ot, then 
expression must always be positive, because the dy- 
namics is described by a divisible map. Therefore for 
any t > and any initial pair of states we must have 
^ > 0. Because the image sets of x{t) and xuit) are 
the same, the condition ^ > must hold also in the 
case of a time dependent decay rate. Stated in another 
way, Eq. ([SJ describes a non-Markovian channel if and 
only if the corresponding dynamical map is non-divisible. 
This condition is valid for any pair of initial states. The 
amount of non-Markovianity of the channel, as defined 
by Eq. ((8)), can be written as 



III. DAMPING MASTER EQUATION 



AAp = max V [.F(P, t+) - tj)] , 



(14) 



We start by considering the dynamics described by the 
following phenomenological Lindblad type equation with 
a single decay channel and a time dependent damping 
rate 7(t), 



dt 2 ^ ' 



a) ap - 



pa^a) 



(9) 



Any Gaussian state evolving according Eq. ([9]) remains 
Gaussian, with the displacement and the covariance ma- 
trix evolving as follows [l3| 



m 

a{t) 



x(t) 

e-^/3(0). 



(10) 



where 



x{t) 



a(0) 



a I 2j{s)ds, 



(11) 



and a being a coupling constant. If x{t) <C 1 then we can 
approximate e"^*-*-* ~ 1 — x(t). Under this weak coupling 
condition we can also approximate the solution of ^ as 
follows 



m 



2 



(12) 



ait)=[l-x{t)]a{0)+x{t)-. 



Upon inserting Eq. (fT2|) into Eq. Q we obtain the ex- 
pression for the fidelity in the weak coupling limit. The 
divisibility property of ([9]) here is equivalent to the con- 
dition j{t) > for any t > 0, as it can be easily verified 
from the solution pUj) . On the other hand, the condi- 
tion for non-Markovianity can be studied by inspecting 



where [t'j',tj] is the I-th negativity interval oij{t). 

For coherent states we may derive an exact expression 
for non-Markovianity valid for any form of the damping 
rate. Assuming a single interval of negativity, we have 

AApe =exp{-A-e-^(*^)}-exp{-Xe-^(*")}, 

\l3ie'^' ~ I32\^ x{t-)~x(t+) (15) 



K = 



-x(t+) 



-x{t-) ' 



where only one angle appears in the parameter K because 
the master equation is invariant under a rotation in phase 
space. Other analytic results can be also obtained if we 
consider small values of the coupling constant. For ex- 
ample in the case of coherent states and for any number 
of negativity periods of the damping rate, we have 



A/pf, ~ max 



Pc 



-x(t+)]K 



maxe^^y^[e-(*?)^-e-(*^")^] 



Pc 



(16) 



a <; max/i(Pc) 

Pc 



2-f{t) dt + o(a^ 



7<0 



where we define the state dependent function /i(Pc) = 
e~^K. From Eq. we can conclude that to first order 
in the coupling a, the quantity A/p^ is proportional to 
the integral of the negativity region of the damping rate, 
and the next contribution appears only to third order. 
The measure is then maximized for if = 1, a condition 
which defines a whole set of pairs of states, all of them 
leading to the same value of Afp^ to first order in a. We 
then conclude that the degree of non-Markovianity for 
the set of coherent states can be written as 



TVpc = -a 7(t) dt + o(a^) 



7<0 



(17) 
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This result has actually been derived under the hypoth- 
esis that also e^'-*^ ^■'^ ~ 1 + x(tf)K for any interval /, 
therefore limiting the parameter domain depending on 
the value of "/{t) and a. However, in the weak coupling 
limit the domain is large enough to contain the points for 
which K = 1. 

We now turn to the class of pure squeezed vacuum 
states (from now on referred to as squeezed states) , char- 
acterized by the set of parameters Ps = {ri,r2, (f>i, (j)2}, 
which, in the absence of displacement, can be reduced to 
Ps = {ri, r2, 0, 0}, with (/) the angle between the squeez- 
ing directions. We can follow the same line of reasoning 
used previously for coherent states and, e.g. for a single 
period of negativity of the damping coefficient [t^,t~], 
define the quantity 



= a \ max5i(Ps) > / 2j{t) dt + o{a) . 



(18) 



where function (7i(Ps) depends on the parameters ri,r2 
and in a rather complicated way. Its expression can be 
simplified if we set r = ri — r2 



(?i(Ps) = 8cosh(2r) 



K^i 0) - \/k{T^ 4>) 



(19) 



where k{r,(j)) = 3 + cos0 -I- cosh4r(l — cost/)). We will 
discuss later on in this Section about the maximization 
of the function 5(1 (Pg). 

From Eq. ([T5)) we notice that the second order term is 
not zero for squeezed states, suggesting that the first or- 
der approximation has a more limited validity than in the 
case of coherent states. It is also nontrivial to determine 
the domain of parameters which allows the expansion to 
be truncated at first order. This is due to the more com- 
plicated functional form of the fidelity when squeezing, 
instead of displacement is implemented. The first order 
result is then of a similar form to the coherent state case, 
with a different state dependent coefficient (7i(Ps). It is 
straightforward to show that this result is independent 
of the number of negativity periods of j(t). 

Interesting results can be also derived if we consider 
now the most general pure Gaussian state (both dis- 
placement and squeezing). From ^ we can separate 
the fidelity as a product of two parts J^{{Pc,Ps},t) — 
C{{Pc,Ps},t)S{Ps,t), where C({Pc,Ps},t) is the ex- 
ponential part, containing the coherent state amplitudes, 
and S{Ps,t) is the fidelity for zero displacement. A first 
order expansion in a shows that 



A/'{Pc,Ps} = a / 27(t)dt 



7<0 

X max {5(Ps,0)/i(Pc) + C({Pc,Ps},0)gi(Ps)} 

{PC:Ps} 



-I- o(a). 



This result gives the non-Markovianity measure at first 
order in a with a pair dependent coefficient given by 
a linear combination of the displacement contributions 
.fi(Pc) and the squeezing one (71 (Ps). However, we can- 
not conclude that even at this order the two contribution 
are completely independent, because the weights appear- 
ing in the combination depend on the fact that we ap- 
plied squeezing and displacement. For example it is easy 
to show that S'(Ps,0) = J^(Ps,0) is the initial fidelity 
of the same pair of states with zero displacement, and 
C({Pc, Ps}, t) = ^({Pc, Ps}, 0)/J-(Ps, 0), is the ratio 
between the initial fidelity and the initial fidelity without 
displacement. 

It is worth noticing that for coherent-thermal states 
with equal thermal parameters, i.e. A^i = N2 — N, at 
first order in the coupling we obtain 



A/'{Pn,Pc} =max 



Pc 



Ar>0 2iV -I- 1 



o(«'). 



(21) 



which is maximized by pure states {N = 0), thus sup- 
porting our choice to restrict the analysis to pure states 
only. 

As we mentioned before many features of non- 
Markovianity for the damping channels do not depend 
on the specific form of the damping rate. However, for 
the sake of concreteness let us now consider an example 
of the damping rate 



1 r e-*/iOsint, 



if 
if 



t < 57r/2, 
t > 57r/2, 



(22) 



(20) 



which is characterized by only one interval of negativity, 
[tt, 2tt], and where we are definitely in the weak coupling 
regime if a < 0.1. 

We start the analysis evaluating the non-Markovianity 
of the channel for coherent states and squeezed states. 
These results are shown in FigdJ where we plot for both 
classes the non-Markovianity evaluated at first order, to- 
gether with the exact numerical solution obtained by 
maximizing ([T4l) using the full solution (fTO|). The co- 
herent states maximization has been carried out over 
all the parameters involved and their domain, while for 
squeezed states we consider two fixed values of the angle 
(j) between the squeezing directions and then maximize 
over the magnitudes ri and r2 , with the maximum value 
reached for ri — r2. Two main results are evident: the 
non-Markovianity of the damping channel is larger for 
squeezed states than for coherent states independently 
of the value of coupling and increases with decreasing 0. 
Due to the nature of the master equation we can also 
state that this result is valid independently of the form 
of 7(i). The second result is that the first order approx- 
imation is well sustained by coherent states while it is 
violated by squeezed states also for relatively small val- 
ues of the coupling. This behavior for squeezed states 
becomes more and more evident if we reduce the value 
of the angle 0, and it can be explained by inspecting Eq. 
for fixed > 0. For small 0, in fact, the function 
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0.02 0.04 0.06 0.08 



FIG. 1: (Color online) Non-Markovianity of the purely damp- 
ing channel as a function of the coupling a. The dotted blue 
line is the result obtained by maximizing over coherent states 
(the exact solution and the first order approximation coin- 
cides). The solid lines are the exact non-Markovianities for 
squezeed states, whereas the dashed lines are the correspond- 
ing first order approximations. Black curves correspond to 
(j> = 0.2,while red ones to (h = 0.1. 



f/i(Ps) decreases, thus determining a violation of the first 
order approximation. A further discussion of the results 
is postponed to the final Section, after we have provided a 
more complete picture by discussing the quantum Brow- 
nian motion model in the next Section. 



IV. QUANTUM BROWNIAN MOTION 

The master equation describing quantum Brownian 
motion (QBM) in the interaction picture, under the weak 
coupling and secular approximations (see, e.g. Ref. [isj 
and references therein), is given by 



dp 
'dt 



+ 



AW±7W.„ t t t ^ 

a (Zapa' — a' ap ~ pa' a) 

A(t) - 7(t) + t t, 

a [2a' pa — aa' p ~ paa ' ) . 



(23) 



The diffusion and damping coefficients A(i) and ^(t) can 
be derived once we provide the analytic form of the spec- 
tral density J{u}) and the temperature T of the environ- 
ment, assumed to be in a thermal state. Their expres- 
sions in the weak coupling limit are 

A{t) — J J rfw J(w) ^A^(w) + i^cos(cjos) cos(ws), 



"y{t) — ds duj J{uj) sm{uJos) sm{ujs), 
Jo Jo 



(24) 



where N{uj) = {exp{huj/kBT}-iy^ is the mean number 
of thermal photons for a mode of frequency w, and wq is 
the bare frequency of the system. The spectral function 



we are going to consider is the following Ohmic spectral 
density 



J{uj) OC UJ C{U!, UJc) , 



(25) 



where C{ui,uic) is a high frequency cutoff function with 
UJc being the cutoff frequency. Usually the cutoff function 
is chosen to be of Lorentzian or exponential form. Here 
we use an exponential cutoff C{uj,uJc) = exp{— oj/wc}, 
such that 



J(a.) 



(26) 



The exact solution of Eq. for the displacement and 
covariance matrix of any initial Gaussian state is 

I3{t) = e-^^(*)/3(0), 

IT /■* (27) 

cr(t) = e-^(*V(0) +a-e-^(*) / e^(")A(s)ds, 
2 Jo 

which, at first order in a can be written as 



2 



m 



(28) 



a{t)^[l-x{t)]a{0) + yit)-, 



with 



x{t)=2a -f{s)ds, y{t)^2a A{s)ds, 
Jo Jo 

and where a{0) and (3{0) are the covariance matrix and 
displacement vector of the initial state, respectively. 

The dynamics generated by Eq. (|23|) is more involved 
compared to that of Eq. and this gives us the op- 
portunity to study in more detail the non-Markovianity 
properties of continuous variable systems. What is here 
relevant, compared to the previous case, is the pres- 
ence of two decay channels, one downward channel and 
one upward. This structure leads to the inequivalence 
of non-divisibility and non-Markovianity. For the mas- 
ter equation (|23p the divisibility property is satisfied if 
A{t) > \^{t)\ [4]. On the other hand, we have a non- 
Markovian behavior if the following quantity attains neg- 
ative values 



dT 



dF dx 
dx dt 



dF dy 

dy dt 



2a-f{t) 



dT 
dx 



-2aA{t) 



dT_ 
dy ■ 



(29) 



The sign of ([29l) depends in a non trivial way on the 
values of the damping and diffusion coefficients, and 
of the derivatives of the fidelity with respect to x(t)^ 
y{t)^ which are in general functions of all the parame- 
ters and the time. Expression (|29|) indicates that both 
damping and diffusion phenomena contribute to the pro- 
cess. The dominance between the two contributions de- 
pends on the spectral function J{uj), the coupling con- 
stant, the temperature of the environment and on the 
initial pair of states through the derivatives dJ- /dx and 
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dF /dy. Each derivative is proportional to the difference 
between two fideUties, the one calculated for an incre- 
ment in one variable and the initial one. We expect that, 
e.g. F[x + h,y) > T{x,y) holds, because, according to 
Eq. p8| , the first term is the fidelity between two states 
which have lost information about their initial prepa- 
ration [see Eq. (1^51) ]. In other words, we expect non- 
Markovianity to be dependent on the sign of the master 
equation coefficients and only on the magnitude of the 
partial derivatives. 

Numerical evaluation of the region of negativity of 
dJ- jdt shows that the zeroes of the derivative are essen- 
tially the same as those of the diffusion coefficient A(<), 
thus suggesting that diffusion is the leading phenomenon 
for non-Markovianity of QBM. In order to prove this re- 
sult we should inspect the form of the non-Markovianity 
at the first order in the coupling a, in the same spirit as 
in Sec. III. For coherent states we have 



e Ja<o 



(30) 



and the contribution from the damping coefficient ap- 
pears only in the third power of the coupling. For 
squeezed states we get 



— = aS'A(ri, r2, 0)A(i) -f- aS^(ri,r2, (j))l{t) + o(a^). 



(31) 



where the coefficients Sj{ri,r2,4') and •S'a (''i , ?'2 j </>) are 
the zeroth order expansions of the derivatives dT/ dx and 
dJ- /dy, respectively. An inspection on the magnitudes of 
the coefficients shows, however, that unless ~ 1 (the 
maximum of the measure is instead obtained in the region 
4> <C 1), we have ^^(ri, r2, (j)) ^ 5'a(?'i, ?'2, 0)- Therefore, 
non-Markovianity can be approximated by 



A/ps = a max S'a {ri ,r2,( 



A<0 



2A{t)dt + o{a'^). (32) 



Within the validity of the first order expansion this ex- 
plains why the zeros of essentially coincides with 
those of A{t). 

As for the damping channel, it is possible to derive a 
closed formula for the non-Markovianity of QBM for co- 
herent states. The only assumption is that the negativity 
of the derivative of the fidelity coincides with the nega- 
tivity of the diffusion coefficient. For a single interval of 
negativity we have 



A/pc = cxp < — 



- exp • 



pg-a;(t+) 



(33) 



where 

2 



3-x(t+) g-^(t ) +y(t-) 



„-x(t-) N -1 



(34) 



-x(i+) e-^(*") -Hy(t-) 



Because the non-Markovianity is here led by the diffusion 
phenomena, we expect a strong dependence also on the 
temperature of the environment. In the following, we will 
examine the non-Markovianity in different temperature 
regimes, comparing the first order approximation and the 
numerical one based on the exact solution (|27l) . 

The temperature dependence of the diffusion coeffi- 
cients is apparent in its definition ([35]) . Essentially the 
coefficient is a sum of a zero temperature term and a 
contribution from the thermal photons in the bath. Re- 
spectively, their expressions are 

1 /■* f°° 
Ao{t) — - ds dujJ{uj)cos{ujQs)cos{uJs), 

^ i° (35) 

Axit) — ds dujJ{Lo)N{uj)cos{uJos)cos{uJs). 
Jo Jo 




FIG. 2: (Color online) Diffusion coefficients as a function of 
time for an Ohmic environment with luc — 1 and (a): cjo ~ 4, 
(b): ujo = 6. We plot the zero temperature diffusion coeffi- 
cient (soUd blue), low temperature {kBT/hwc = 1/5, red dot- 
ted), intermediate temperature {kBT/huJc = 1, black dashed) 
and high temperature {kBT/hwc — 4, green dotted-dashed 
line). 

For low temperature, that is when the thermal en- 
ergy k^T is much smaller than any excitation energy, 
hujQ, hide, the total coefficient does not differ much from 
the zero temperature contribution. Therefore we expect 
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non-Markovianity to be independent of T in this regime. 
As T is increased ^xit) starts to become relevant and 
in the high temperature regime becomes dominant com- 
pared to Ao(t). In this regime the coefficient depends 
hnearly on T and therefore we can expect a strong de- 
pendence of non-Markovianity on temperature. Another 
feature is ihustrated in Fig. [21 where we show the diffu- 
sion coefficients for different values of the bare frequency 
wq. By comparing the two panels, we can see that if 
we approach the resonance condition Wc ~ Wq, the dif- 
fusion coefficient at low temperature does not show neg- 
ative regions and so the first order non-Markovianity is 
vanishing. As we increase the temperature, A(t) may 
become negative and therefore the system shows a non- 
Markovian behavior. This distinction between low and 
high temperature ceases to be valid when we are out of 
resonance ujc/uJa <C 1, where also for low temperature we 
can have a non-Markovian behavior. 




FIG. 3: (Color online) Non Markovianity of QBM channel 
with Ohmic reservoir on coherent states as a function of the 
coupling a with ujc ~ 0.2 and ujq ~ 1, for different tem- 
peratures. Dashed (solid) lines represent the first order ap- 
proximate (numerical) solutions. Black curves correspond to 
fcflT/^o = 0.2 and blue curves to kBT/hwo = 0.5. The inset 
shows the temperature dependence for fixed coupling a = 0.1, 
and cjo = 1 (red solid), cjq = 1-1 (red dashed). 

After this discussion on the nature and behavior of the 
diffusion coefficient, we are ready to illustrate the results 
about non-Markovianity of the QBM channel. We start 
by considering the class of coherent states, for which the 
first order measure is defined in Eq. ([5(11) . In Fig. [31 we 
show the comparison between the analytic and numeri- 
cal results for ujc = 0.1 and ljq — \, and for two different 
values of temperature. When a < 0.1 the agreement is 
good in both cases, whereas the first order approxima- 
tion fails for the higher temperature when we increase 
the coupling. In the inset of Fig. [3l we show the first or- 
der measure for fixed coupling (a = 0.1) as a function of 
the temperature when Uc — 0.3 and different values of 
ujQ. When T < 0.4 the measure is zero indicating that 
the diffusion coefficient has no negativity regions. 

For squeezed states, in Fig. [31 we show the comparison 
between the numerical and first order analytic results for 



the measure as a function of the coupling constant a for 
small values of the temperature. We can notice a behav- 
ior similar to that of the damping channel case (See Fig. 
2a), the exact and approximated expressions indeed co- 
incide only for very small value of the coupling. 




FIG. 4: (Color online) Non-Markovianity of QBM chan- 
nel with Ohmic reservoir as a function of a. We have set 
ksT /fvjjQ = 0.2, = 0.2, uiQ = 1. Red and black curves show 
the behavior of squeezed states, blue curves are for coherent 
states. Dashed (solid) curves represent the first order (nu- 
merical) solution to the maximization problem. For squeezed 
states the maximization has been performed with the con- 
straint ri = r2 and for fixed phase </> = 0.1 (red curves), or 
cj> — 0.05 (black curves). For coherent states the maximization 
has been performed over all their characterizing parameters. 

For increasing a the non-Markovianity saturates to a con- 
stant value, which is achieved for smaller a when (p de- 
creases. In Fig. Ulwe compare the results for coherent and 
squeezed states, showing that also for the QBM channel 
squeezed states are more sensitive than coherent states 
to the non-Markovianity of the channel. 




FIG. 5: (Color online) Non Markovianity of QBM channel 
with Ohmic reservoir on coherent states. We have set cjc ~ 
0.2, ujQ — 1 and maximized at fixed angle (p = 0.05 over the 
squeezing parameters. Different lines correspond to different 
temperatures: kBT/hwQ — 0.3 (black solid), kBT/hwo = 0.9 
(blue dashed), kBT/huiQ = 4 (red dotted) and kBT/huio — 8 
(green dotted dashed). 
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Finally in Fig. [5] we plot for fixed parameters Wc = 0.2, 
= 1 and (j) — 0-05 the numerical results for the mea- 
sure for different temperatures. The behavior is quali- 
tatively the same, i.e. we have saturation for increasing 
a and the saturation value increases with temperature, 
until it reaches a maximum saturation value for high tem- 
peratures. 

V. DISCUSSION AND CONCLUSIONS 

In the previous Sections we have analyzed in detail 
the non-Markovianity of two kinds of CV quantum chan- 
nels, addressing separately the non-Markovian behavior 
for coherent and squeezed states. The dynamics of co- 
herent states is governed by the evolution of only the 
displacement amplitude in the damping channels, and 
by both the displacement and the covariance matrix in 
the QBM channel case. On the the other hand, the dy- 
namics of (zero amplitude) squeezed states depends on 
the covariance matrix only in both cases. This behavior 
allows us to address the effects of the two contributions 
on the amount of non-Markovianity of the channel. 

Despite the difference in the form of the master equa- 
tions, many common features are apparent. At first, we 
notice that squeezed states are in general more sensi- 
tive to the non-Markovian behavior, as witnessed by the 
larger non-Markovianity compared to coherent states, in- 
dependently on the value of the coupling and temper- 
ature. Generally, for both classes of states, the non- 
Markovianity is mostly related to the sign of the mas- 
ter equation coefficients. For the case of QBM non- 
Markovianity is mostly due to diffusion described by the 
coefficient A(<). On the other hand, even if the dynam- 
ics of the displacement vector is not affected by it [see 
Eq. p7)) ]. it is still fundamental for the class of coherent 
states. This behavior is even more evident as the tem- 
perature is increased, since the damping is independent 
of temperature. 



Another interesting feature is the behavior of the non- 
Markovianity for squeezed states as a function of the cou- 
pling and the temperature, in particular for what con- 
cerns the quantum Brownian motion case. As it is ap- 
parent from Fig. [5l the value of non-Markovianity as 
a function of the coupling saturates, with a saturation 
value that increases with the temperature. This behav- 
ior, whose origin may be traced back to the time evolu- 
tion of the covariance matrices, implies the existence of 
some bound on the flow of information from the environ- 
ment back to the open system as a result of the Gaussian 
structure of the map. 

In conclusion, we have introduced a measure to quan- 
tify the non-Markovianity of continuous variable quan- 
tum channels and have used it to analyze two paradig- 
matic Gaussian channels: the purely damping channel 
and the quantum Brownian motion channel with Ohmic 
environment. We have considered different classes of 
Gaussian states and found the pairs of states maximizing 
the backflow of information. For coherent states we have 
found analytical solutions, whereas for squeezed states 
we have resorted to numerical maximization, and also 
obtained some approximate analytical solutions in the 
weak coupling limit. 

Our results are encouraging enough to suggest the use 
of our measure of non-Markovianity to analyze more gen- 
eral Gaussian channels, and to assess non-Markovianity 
as a resource for quantum technologies. 
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